clear all
clc
ft = @(x)sin(x)/x;
I= mymulNewtonCotes(ft,0.01,5,10,2)